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Abstract 


We propose a novel formalism for simultaneously describing both, the hard 
and soft parton dynamics in ultrarelativistic collisions of nuclei. The emission 
of gluons from the initially coherent parton configurations of the colliding nu¬ 
clei and \ow-pt color coherence effects are treated in the framework of a Yang- 
Mills transport equation on a coupled lattice-particle system. A collision term 
is added to the transport equation to account for the remaining intermediate 
and high-pt interactions in an infrared finite manner. 


Experiments with heavy ion collisions at energies above 100 GeV/u, in preparation at 
the Relativistic Heavy Ion Collider (RHIC) in Brookhaven and the Large Hadron Collider 
(LHC) at CERN, will try to establish the existence of a new phase of nuclear matter, the 
quark-gluon plasma |^. One of the theoretical challenges in this context is to develop a 
description on the basis of quantum chromodynamics (QCD) of the processes that may lead 
to the formation of locally equilibrated superdense matter in these nuclear reactions. 

In recent years, most theoretical attempts at developing such a description have been 
based on the idea that, at very high energy and for heavy nuclei, the dominant mechanism of 
energy deposition in the central kinematical region is the perturbative scattering of partons 
M- Because the interactions among gluons are stronger than those involving quarks, 
this mechanism predicts an abundance of gluons during the early equilibration phase 0. 
This concept can be generalized into a theoretical framework, called the parton cascade 
model 1^, which formulates thermalization as a transport process involving perturbative 
QCD excitations, i.e. quarks and gluons PJ. The predictions of this formalism have been 
extensively studied by means of numerical simulations [MITU . 


One of the problems inherent in this formulation concerns the description of the initial 
state. The transport equations start with the assumption of a probabilistic phase space 
distribution of partons, whereas in reality the states of the colliding nuclei are described by 
coherent parton wave functions. The incoherent parton description fails, especially, at small 
transverse momenta, because the QCD coupling constant diverges in naive perturbation 
theory. Some time ago it was proposed that the proper solution to these difficulties would 
be the perturbative expansion, not around the ’’empty” QCD vacuum, but around a mean 
color held describing the static color held accompanying the fast moving valence quarks of 
the colliding nuclei [|^ . 
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Because the mean color field of a heavy nucleus locally receives contributions from the 
quarks contained in many different nucleons, its source can be represented as a Gaussian 
ensemble of color charges moving along the light cone |^. We will, therefore, refer to 
this model here as the random light-cone source model (RLSM). Within this framework, 
the energy deposition by gluonic interactions is described as classical gluon radiation at 
small transverse momenta and as gluon-gluon scattering at high transverse momenta 

r5| . Quantum corrections to this picture predict an enhancement of the glue field of 
the colliding nuclei at small values of the Bjorken variable x. The full solution of the 
nonlinear classical RLSM equations for the color field of two colliding nuclei requires a 
lattice formulation in 2-1-1 dimensions |T7[ . 

The possibility of a description of inelastic gluon processes by means of the nonlinear 
interactions of classical color fields has also been explored numerically in studies of collision of 
two Yang-Mills field wave packets on a one dimensional gauge lattice . These calculations 
gave evidence that the interaction between localized classical gauge fields can lead to the 
excitation of long wavelength modes in the collision, which is reminiscent of the production 
of an equilibrated gluon plasma. 

Here we address the question how this new insight can be incorporated into the con¬ 
ceptual framework of the parton cascade model. First of all, it is necessary to include a 
coherent color field in addition to the incoherent quark and gluon distributions, qi{r,p) 
and g{r,p). The subscript “f’ here denotes the various quark flavors. We will also insist 
on a full (3-|-l)-dimensional representation, which will permit the study of deviations from 
boost invariance. 

Because even the classical Yang-Mills equations do not, in general, allow for global an¬ 
alytic solutions [T^J2^, we propose to solve the RLSM equations numerically on a gauge 


lattice. Lattice calculations in Euclidean space-time have been shown to provide a reliable 
approach for the calculation of static and quasi-static properties of strongly coupled quan¬ 
tum field theory, in particular, QCD. For dynamical systems far off equilibrium, however, 
one needs to study the system in real continuous time. The lattice discretization then should 
only be applied to the Euclidean sub-space R^. In this case it is appropriate to choose a 
Hamiltonian formulation rather than a Lagrangian one. We have to emphasize that this 
concept is neither explicitly invariant under general gauge transformations nor Lorentz in¬ 
variant. However, we believe that for the type of problems described above this method is 
indeed useful. 

One has to select a rest frame in the space R ® R^ which in our case probably is best 
chosen as the center of velocity. Further, one has to adopt a gauge. The temporal gauge 
in the continuum = 0) seems most appropriate here |^. A set of equations describing 
the evolution of the phase space distribution of quarks and gluons in the presence of a mean 
color field, but in the absence of collisions, was proposed more than a decade ago by Heinz 
2^,|23[]. This non-Abelian generalization of the Vlasov equation can be considered as the 


continuum version of the dynamics of an ensemble of classical point particles endowed with 
color charge and interacting with a mean color field. The equations for this dynamical 
system were originally derived by Wong ||2^ . 


In the following we develop a formulation of the RLSM including the ideas of Heinz and 
Wong. We represent the valence quarks of the two colliding nuclei as point particles moving 
in the space-time continuum, and interacting with a classical gauge field defined on a spatial 
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lattice but with quasi-continuous time.[] In principle, this idea follows the proposal of Hu 
and Miiller [^| for the simulation of the effects of hard thermal loops by means of colored 
point particles. 

At this stage we are still general enough to assume that the soft modes of the gluon 
helds are described through gauge helds with SU(iV) symmetry. In the associated Lie- 
algebra LSU {N) we express the Hamiltonian of the above outlined system in the continuum 
as 


Ni _ N2 _ 

^ = J2\/\Pi\^ + ^0 + Y. \/\Pi? + ^0 

1=1 i=Ni+l 

-2gj £xi:i[J^A^] - (1) 

where A^, Qt e LSU(A^). The curly quantities denote those in the adjoint represen¬ 
tation, which are defined as, e.g., A^ = ■ T'^ with group generators T'^. g is the gauge 

coupling constant, and denotes the held strength tensor of the mean color held The 
moving particles generate a color current 

= 5 A 3 -mi ( 2 ) 

i \J\Pir + m^ 

where t denotes the global time in the chosen reference frame. Denoting the space-time 
positions, momenta, and color charges of the particles by xf, pf and respectively, the 
following equations of motion are derived from the above Hamiltonian (|^): 


odXi ^ 

Pi-r- = Pi 
dt 

( 3 ) 

(jr)^ 

p<^^ = 2gTT{Q,J^nP^,u 

( 4 ) 

p^i Pfx,i 

( 5 ) 


The factors on the l.h.s. are needed to convert the derivatives with respect to proper time 
into coordinate time derivatives. Furthermore, the inhomogeneous Yang-Mills equations 


D,T>^^{x) = gj'^ix) 


( 6 ) 


describe the dynamics of the classical mean color helds. The current density (^) forms 
the source term on the r.h.s. of (^. The coupled system of the Wong equations (M 
and the Yang-Mills equation (^) is highly nonlinear and can only be solved numerically or 
perturbatively. 

These equations have been used to simulate the ehects of hard thermal loops on the 
dynamics of soft modes of a non-Abelian SU(2) gauge held at hnite temperature p7|,|^. In 


^The numerical implementation also requires a discretization of the time variable, but the temporal 
step size can be taken arbitrarily small. 
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this case, the colored particles describe the gauge field modes with thermal momenta, and 
the mean field describes the coherent motion of those gauge field modes which have a wave 
number k much smaller than the temperature T and are highly occupied. The separation 
of the two regimes was achieved by discretizing the mean gauge field on a lattice with 
elementary spacing a -C T~^. Requiring particles to have momenta p > ir/a then avoids 
double counting degrees of freedom. 

Here we propose to use the equations (@-§P to describe the interactions among the glue 
field components of two colliding heavy nuclei. In this case, the lattice cut-off a can be used 
to separate the regime in transverse momentum where the dynamics of gluons is perturbative 
(large k^) from that where naive perturbation theory fails (small /ct)- The gluon propagators 
used for the calculation of the collision terms will be regulated in the infrared by the lattice 
cut-off kc = 71/a. The interaction with the mean color field allows for an exchange of an 
arbitrary number of gluons, and the screening of the soft components of the gauge field by 
perturbative partons is taken into account naturally by the nonlinear nature of the 

coupled equations (M). 

Following the idea of Kogut 


we approximate the gluonic part of the Hamiltonian 


(0) by a discretized form on a gauge lattice. In contrast to Ref. |^T|, however, we represent 


the fermions through point-like particles. This leads to a Hamiltonian which is represented 
as a sum of the following terms 


H= H, 


part 




(lattice) 

YM 


(7) 


where impart contains the first two terms on the r.h.s. of (|l|) and is defined as 




(lattice) 


YM 


= —a 


^ Tr| Sx,kSx,k 


x.k 


^ l,m 

9 J^x,k 


( 8 ) 


As already mentioned, the dynamical equations (^|-^ can be solved efficiently by nu¬ 
merical time integration. A lattice version of the continuum equations is constructed |]^J^ 


by expressing the gauge fields in terms of link variables Ux,ie SU(A^), which represent the 
parallel transport of a field amplitude from a site x to a neighboring site (x -f /) in the 
direction /. As in the Kogut-Susskind model 


2^ we choose the temporal gauge Aq = 0 and 


define the following variables: 


Ux,i = exp{-igaAi{x)) = 

Ux^kl Ux^kUx-\-k^l^x-\-k-\-l^—kUx-\-l^—l 


(9) 

( 10 ) 


Consequently, we have 


P ■ — ^ 77 /yt 

~ iga 


Aiga^ 


^jkl {^x,kl ^x,kl) 


( 11 ) 

( 12 ) 
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for the electric and magnetic fields {Sxj, Bxj e LSU(iV)), respectively. There are advantages 
in choosing Ux^i and £x,i as the basic dynamic held variables. This choice transforms the 
discretized Yang-Mills equations into the following equations of motion 


4<=(*) = - K.mV) (14) 

- uiumUkiim-iAt) 

9 £^x,k(t)- 


In the spirit of the statistical nature of the transport theory, we split each quark into a 
number Ug of test particles, each of which carries the fraction qq = Qojug of the quark color 
charge Qq. In a hrst step, we adopt the gauge group SU(2) here for simplicity. Consequently, 
each nucleon is represented by two quarks (instead of three), initially carrying opposite color 
charge. 

Perturbative short range interactions at high momenta can be described in form of a 
stochastic collision term, well known from Boltzmann-type transport equations [pl|,^. For 
a consistent description of both, long range and short range interactions on an equal footing, 
the equations of motion ^ for the long range interactions have to be cast into the form 
of a single transport equation and combined with the collision term. The Vlasov part of the 
transport equation, from which the equations of motion (|-i can be recovered, was hrst 
derived in We extend the formulation by adding a stochastic collision term similar 

to the one used in p. The full transport equation then follows as: 


+2tgTT{[Qi,A^]-dQ)}ik{Xi,Pi, Qi) 

= Y. (15) 

processes 


Here f^ denotes the one-particle distribution functions of the valence quarks and of the 
“hard” gluons {k = q, g). This set of nonlinear integro-differential equations is coupled 
to the Yang-Mills equation in which the color current is now given by a moment of the 
one-particle distribution functions: 


= s E /da % aprt,(pr,<, a) 

k '' Pi 

The collision integrals have the form: 

(^(pr.xf.a.r) = ^ ■/«(|p.l - l*:ol) ndF;/ 

{2n)‘ D{h(tA, 4. a)) 


(16) 


(17) 


with 


5 








(18) 


D(f»(pr.<.a)) = JlWuxt.Qi) - nfs(p.‘'.<.a 


out 


and 


TT dr ■ = IT __ 

^ (2vr3)(2p0) 


^(bil - \kc\) • 


Si is a statistical factor defined as 


(19) 


Si = n K'. KT'. (20) 

j¥=i 

with identical partons of species a in the initial or hnal state of the process, exclnding 

the ith parton. 

The step fnnctions 9{\pi\ — \kc\) insnre that only hard particles are allowed to propagate 
in the system. The snperscript (c) on the matrix element Ai indicates that only the hard, 
i.e. short range, part of the interaction is treated in the collision term. This cnt-off will be 
discnssed in more detail below. 

The matrix elements acconnt for the following processes: 


A q + q'^ q + q' , 
B q + q^ q + q , 
C q + q^ g + g , 
D g + g^ g + g , 


( 21 ) 


together with those obtained from crossing relations {q and q' denote different qnark flavors). 
The amplitndes for these processes - not taking the infrared lattice cnt-off kc into acconnt 
- have been calcnlated in refs. [p3H35[| for massless qnarks and in refs. for massive 

qnarks. The corresponding scattering cross sections are expressed in terms of spin- and 


color-averaged amplitndes 

di 


Ibvrs^ 


(|Al<“>(J,t,fi)R 


( 22 ) 


with s, t, u being the well-known Mandelstam variables. For the transport calcnlation we 
also need the total cross section as a fnnction of s which can be obtained from (^): 



(23) 


The integration bonndaries are hxed throngh kinematical constraints. Note that the treat¬ 
ment of the cross section (pT|-[23|) with the matrix elements snpplied in [p3|-|^ does not take 
the infrared lattice cnt-off kc into acconnt. The rigorons way to evalnate the matrix elements 
lAf P and to eliminate the small momenta from the glnon propagators wonld be to snbtract 
the lattice propagator from the continnnm propagator in the Feynman diagram describing 
the scattering process at lowest order. Becanse the evalnation of the glnon propagator on 
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the lattice is complicated, we propose here to use, for exploratory studies, the usual matrix 


elements [^^ 37| but with a cut-off on the allowed momentum transfer, corresponding to the 


lattice cut-off kc = Tr/a. We can cast this into the Lorentz invariant form that the scale of 
the interaction, Q‘^{s,t,u) must satisfy the constraint 


Q‘^{s,i,u) > kl 


(24) 


The functional form of is generally process dependent and not unambiguous, although at 


leading order all choices for Q that increase with the parton-parton center-of-mass energy are 


equivalent. One can now solve equation (|2^) for t in order to obtain an additional constraint 



for the integration boundaries of equation (^). Thus, only momentum transfers larger than 

kc contribute to the total cross section. It was shown in ref. that spectrum of the classical 
Yang-Mills radiation matches smoothly onto the conventional minijet distribution near the 
intrinsic transverse momentum scale of the partons in a heavy nucleus at high energy. The 
resulting expectation that the precise choice of the momentum cutoff kc is not important 
must, of course, be verified by future numerical calculations. 

In summary, we have developed a novel formalism, which allows for the first time the 
treatment of both, the hard and the soft parton dynamics in ultrarelativistic heavy ion 
collisions in a consistent transport approach: The emission of gluons from the initially 
coherent parton conhgurations of the colliding nuclei as well as low-pt color coherence effects 
in parton-parton scatterings are treated in the framework of a Yang-Mills transport equation 
on a coupled lattice-particle system. Intermediate and high-pj’ interactions are described in 
a collision term similar to that of the parton cascade model. This formalism thus avoids 
problems connected to the infrared cut-offs in the parton cascade model and offers a unihed 
treatment of coherence effects within that approach. 
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